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применяемых в сельском хозяйстве 


Farm organization planning calls for taking into account many 
variants of possible solutions of the problem as well as for adjustment 
to many constraints imposed by natural and economic conditions. Due to 
this fact, linear programming is recognized as an efficient instrument of 
optimizing production and investment plans, although not free from 
defects. One of them is that conventional linear programming cannot 
adequately cope with fluctuations of crop yields, prices and of other 
parameters. Consequently, there was a long-prevailing opinion that 
linear programming could be used in very rare cases only. Over the last 
several years, however, significant progress has been made in the 
so-called stochastic programming, especially in its theory. Numerous 
methods also appeared which could be, and indeed were, applied in agri- 
culture. It would be useful to present at least some of the most impor- 
tant formulations. For the interest in linear programming is relatively 
high while there is little information in Polish scientific literature on 
the methods of risk considering. 


1. FORMULATION OF THE PROBLEM 


A standard version of the linear programming problem is the follow- 
ing: 

maximize mTx, such that: 

Ax<b and x>0 
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where: 

mT — a column vector of objective function parameter mean values, 
x — vector of activities, 
A — an input-output coefficients matrix, 
b — vector of available amounts of scarce resources. 


We have to assume that the vector m and also the matrix A are 
subject to fluctuation because fluctuations of crop yields and prices 
cannot be excluded. In some cases the vector b has to be considered as 
well: the amount of available labour in respective periods of different 
years can differ due to changing weather conditions. 


II. RISK CONSIDERING IN OBJECTIVE FUNCTION PARAMETERS 


If we assume that prices are the only source of fluctuations of agri- 
cultural planning parameters or that all crops grown cn the farm are 
cash crops, it is sufficient to concentrate on objective function parame- 
ters only. This assumption, apparently artificial in farm conditions, is 
useful to the extent that it permits to see the approach to the problem 
of risk in objective function parameters. The extension of chance action 
upon other elements of the linear programming model, that is an input- 
-output coefficient matrix and a right-hand side vector, does not in any 
way affect the approach to the introduction of risk into the objective 
function. In the two oldest and best-known methods of considering 
risk in objective function parameters, formulated by Markowitz (14) and 
Freund (10), the measure of fluctuation is the total variance of objec- 
tive function: 


Var, = XTDx,! where: 
D — variance-covariance matrix of objective function para- 
meters, 
xT, x — column and row vectors of activities, respectively, 
Vm x — total variation of objective function. 


The Markowitz method (14) was originally meant for choosing 


. n n 
1 Or another way: Ve = Z xo; + У Exxon where 
j=1i=1 
xi — i-th activity, 
xj — j-th activity, 
с — variance of objective function parameters of the i-th activity, 
бі — covariance between parameter values of the i-th and j-th activities, 
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a stocks combination, hence its name of "portfolio selection”. It is 
founded on the assumption that the goal of financial activity is to maxi- 
mize the profit, which, translated into formulas of linear programming, 
means a maximization of mTx objective function, with Ax « b and 
x > 0. The Ax <b constraints are necessary because the amount of 
money that could be spent for stocks by any individual or company is 
limited just as is the amount of a single firm's stocks available on the 
market. It also follows from this assumption that the maximization 
should: be such that the total variation of profit dces not exceed a cer- 
tain value, which could be accepted by a decision-maker. This means 
that an additional constraint has to be imposed on mTx, that is xTDx < 
< a, where a the maximum admissible value of objective function varia- 
tion. Since this is an entirely subjective value and it is difficult to assu- 
me any relation between a and mTx in advance, the most convenient way 
of solving this problem is to use parametric programming, with the pro- 
blem formulated as follows: 
maximize: mtx, such that 

Ax $b 

xTDx За 

x 20, where: 3 


тт, x, A, b, D and a as above. 


Such problems could not be solved in the early fifties when the 
”portfolio selection” method was formulated. The converse problem, that 
of minimization of objective function variation, with the assumption that 
the mean value of the profit will not decrease below a certain value, 
could already be solved owing to the earlier work by Kuhn and Tucker 
(12). Its mathematical solution turned out to be identical with the original 
problem. The final version of the Markowitz method can thus be for- 
‘mulated as follows: 


minimize: xTDx, such that: 
Ах < b, 
m™x > В, 
х 20, where 


m, x, A, b, D as above 
ß — parameter determining the minimum acceptable 
| profit. : 


The solution to this problem are pairs of mean profit value and profit 
variations, and a corresponding set of values of each activity involved. 
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A choice is made according to individual preferences of profit height 
and its variance. In other words, solutions after the "portfolio selection” 
method provide information that with a given mean profit value, varia- 
tion equal to xTDx cannot be avoided, and that in this case all activities 
have to asume the values as in the optimum solution to the foregoing 
problem. . 

Freund founded his method on the "utility theory" formulated by 
von Neuman and Morgenstern (9). The central point of this theory is 
the assertion of a decrease in money value following its acquisition 
uncertainty. This means that of two farm enterprises with the same 
amount of profit, the one with a lower profit variation is "more useful". 
Moreover, two enterprises with different profits and with a different 
rate of profit variation have equal "utility" if the enterprise with 
a higher profit variation obtains this profit higher by a definite amount. 
This vàlue varies with every individual farm operator. The relative 
measure of this value is referred to as a "risk aversion coefficient". The 
relation between profit height and its variation and profit utility is 
called utility function.” The one proposed by Freund for farmers has 
the following form: 


f(u) = xe PH, where 


e — natural logarithm base, 
a — risk aversion coefficient, 
r — profit height. 


The bigger a is, the less readily a farm operator will take up risk, 
end the higher profit has to be obtained to level higher variation in 
alternative activities. Assuming r to have a normal distribution, the 
expected utility value will be as follows: 


E(u) = p—ao*/2, where: 


u — mean profit value, 
a — risk aversion coefficient, 
o? — profit variation. 


Translated into linear programming, this means: 
maximize: mTx— > xTDx such that: 
Ax <b 


x20 
4 


2 Freund called this relation the "utility of money function”. In other papers, 
the term "utility of function” can be encountered. 


Methods of Considering Risk in Programming Models... 367 


Solution of the above problems, where risk has been dealt with 
according to the two presented methods, requires quadratic program- 
ming. Available computers solve that problem easily, nevertheless, 
quadratic programming is far less convenient than linear programming, 
mainly because the size of problems'is then much more limited. Hence 
there were attempts to modify and adapt the portfolio selection and 
Freund method to the simplex procedure, and to linearize the objective 
function. | - 

The best-known linearization of the portfolio selection method is 
the so-called MOTAD proposed by Hazell (11). Its guiding idsa is to 
replace variation by absolute deviation. Hazell assumes further that it 
is sufficient to take into account negative deviations only. The resulting 
formula is as follows: 3 


minimize Y yr such that: 


ігі 


. 
Ax <b 

m'x > В, where 

yy — negative absolute deviation of j-th activity 


from its mean profit value. 


Chen and Backer (7) prcposed a linearization of objective function 


a 
E(u) = m"x- x'Dx, founded on the assumption that no activity can 
be activated beyond the point where its marginal utility assumes a zero 
value. This means that that the value of any activity can be increased 
as long as its increase adds anything to the sum of total utility. If this 
limit is exceeded, total utility decreases. This marginal utility equals: 


à E(u) 
ÒX; 


a n 
= m— — x S; X;, Where: 
25 


s; — covariance between the i-th and j-th activities 3. 


3 The value of convariance between the objective function parameters of the 
i-th and i-th activity is the variance of objective function parameter of the i-th acti- 
vity, with the marginal utility thus being as follows: 


oE(u DA e : | 
w = m, = aS;x,-L.2a X $нк» where 5; із of course the variance 


j=1 


OXI 


of objective function parameter of the i-th activity. 
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Consequently, the problem to be solved is as follows: 
maximize: m'x, such that 


Ax <b 


n 
4 


a 
шараға "= 

ізі UON 
wherei=1..n,j=l..n 


n x the number of activities 
: ir) 


Or another way: g^ Se 


maximize: m'x, such that: 


whoa 
Ax <b v Top i 
2 W 
Dx < —m u" Ы 
а "ox 
x 20 tolse CIE E 


um 
„i 


Unfortunateiy, so simply formulated a proi ‘gm can be solved only if x 
assumes a positive values. If any activity which, the x vector consists, 


say Xy, assumes a zero value, it could tu.n ` that the constraint: 
а qu 
my — > >> 5% Х; > 0 


ізі 
is restrictive to other activities although x, should” not “tave any 
influence on the optimum solution, because it is an idle activity. Chen 
and Backer developed a multi-stage algorithm for this purpose, which 
gradually removes all idle activities and their corresponding con- 
straints which ensure the assumption of its non-negative marginal uti- 
lity. The algorithm is the following: 


1. Find an optimum solution of a parametric L.P. problem: 
maximize: m'x, subject to: 
Ax <b 


2 
Dx < — m 
а 


x >20 where: 


a — a parameter assuming values from +оо to 0. 
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2. Record the whole set of solutions and their objective function 
values, if none of the dual solutions associated with the constraints 
which are to preserve non-negative x; utility, assumes a positive value. 

3. Remove from the x vector all the activities which are not in the 
basis and all the corresponding constraints ensuring non-negative x, 
utility. 

4, Find a new set of solutions. Come back to step 2 and record only 
ov.dions with a lower m?x valué than previously obtained. 

Another way of the linearization of objective function in the Marko- 
«tz method is "Separable Programming" (27). 

It consists in the division if the xTDx function into a sum of single- 
argument functions, which pei. .ts their spatial linearization. 

The third approach to the pvo of objective function fluctuations 
is founded on the theory of game, . the Polish economic literature 
44s approach has been described i: le b; T. Marszalkowiez (15). 
It appears, however, that a further d.:cussion will be more lucid if the 
basic tenets of the theory are explained з 15 point. | 

In the farm organization planninr :« other decision making, the 
sc-called games with na. u. are s мей out of a number of games 
covered by this theory. 1.- se. Е’ ies have such a property that the 


opponent in the game — na ur: — although ruthless, is not spiteful. It 
is therefore assumed that a ".yer — in this case a decision maker — 
faces m possibilities, each Фет having n realizations of the value 
under consideration. Th: :m is to select one out of m possibilities, 
the choice in no way the opponent’s action. The selection need 


not be limited to the c. „of one possibility, which is called "pure stra- 
tegy". This can also be;any combination of possibilities, which is then 
called "mixed strategy"i Mean values or variances as a criterion of 
choice cannot be applied; as they cannot be calculated because nothing 
is known about the probability of any m realizations.* The only informa- 
tion we have is the set of values which every m possibility can assume. 
To deal with this really difficult situation, the minimax rule is adopted 
if the realizations of m are costs, and the maximin rule if the realia- 
tions of m are incomes. 

The minimax rule consists in the choice of such a pure or mixed stra- 
tegy that has the lowest maximum cost value of all mixed and pure 
strategies. By analogy, the maximin rule selects such a mixed or pure 


4 The mean value as a criterion of selection is called the Laplace criterion. It 
is based on the assumption that if the probability of no m realizations can be de- 
termined, it is necessary to assume that the probability of each realization is the 
same. This approach has been criticized in paper (24). 
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strategy that has the highest minimum value of income of all mixed 
or pure strategies. 

An example will serve as a better illustration. In Table 1 are shown 
the yields of four oats varieties in the course of five years. The results 
Are given in pounds per acre (the example was drawn from Heady, 

Pesek and Walker (29). l 

According to the maximin rule, the B variety is the best because its 
lowest yield obtained in the first year $f the experiment is higher than 
the lowest yield of any other variety under consideration. 

Choosing a mixed strategy is much more complex. For it is impossible 
to make a set of all combinations since the number of proportions of 
each variety in such a combination is infinite, whereas the mixed sta- 
tegy is supposed to have such proportions of each variety that a combi- 
nation with a higher minimum yield could not be found. It is therefore 
necessary to solve the following L.P. problem: 


maximize: xs, such that: 


1472x; + 1568x2+ 1440x; +1552x4— x, > 0 
2112x,+1984x2+2368x3 + 2688x4—x5 > 0 
1920x,+1824x>+ 2496x3 + 2784x4—x5 > 0 
3620х, 1 3104х»- 3552x3+  Oxy—x; 20 
3072x, +3328x, + 2848x3 + 3200x4—X5 > 0 


xit xt x3+ Ха =1 
Хі > 0 
Хә > 0 
Xg 20 
X4 > 0 
Х5 > 0 


The solution to the above problem is a mixed strategy consisting of 56 
of B variety and 44°/o of C variety. 

The whole problem of determining a mixed strategy can be genera- 
lized as follows: 


+ Xp maximization, such that 
M¡X—Xp > 0 


mx— Xp > 0 
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m,Xx— хр > 0 
DIE <d 
isl 
x>0, where: 
Xp — value of a game, 
m; ... M, — vectors of m realizations, 
x — activity vector, 
d — value which the sum total of activity values cannot 
exceed (most often 1 or 1009/o), 
x; — constituent activities of the x vector. 


After this theoretical discussion, it is necessary to return to the appli- 
cation of the method in the construction of an LP matrix which is to 
determine the optimum program of production and possible investments. 
The set of constraints constituting the mixed strategy contains an 


n 
DIE <d element. In the matrix constructed for the described task, the 
=] 1 
> x xd is replaced by the whole input-output coefficients matrix. The 
ici Я 
problem can thus be formulated as follows: 

maximize: Xp, such that: 

Ax <b 


m,x— x, > 0 


The criterion of choice used in the foregoing example is not the only 
one, although the most popular. A detailed analysis of all criteria can be 
found in Adamus (1). 

A similar approach to objective function fluctuations as in the theory 
of games can be found in the "safety-first" method. The idea of "safe- 
ty first” was worked out by Roy (26) and Tesler (28). It was further 
developed and applied to LP by Maruyama (16) and by Petit and 


24* 
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Boussard (21). According to this method, a farm should be operated 
in such a way that the profit every year could be high enough for 
the farm to maintain its existence. This means that the farm’s income 
has to ensure at least a social minimum for the farmer and his family, 
and to pay for all the charges (debet installment payments, interests, 
taxes etc.) every year irrespective of weather conditions and price 
fluctuations. It is not enough to have a high mean income because it can 
be spent if "bad harvest” is not expected. Moreover, previous incomes 
do not necessarily imply that they will be similar in the future. At best, 
it only follows that such and such incomes, profits, or yields will be 
obtained in the future. It is impossible to know how often this will 
happen for the sample is too small to infer anything from, the more 
so that the observations from the previous years cannot possibly be re- 
cognized as drawn out by lot. 

Therefore, the LP matrix should be such as to prevent a situation 
where the mean profit or income is high, but its stability is not suffi- 
cient, which leads to a farm failure. In Maruyma’s already-cited work, 
this problem is solved by: 


mx maximization, such that: 
Ax <b 


mx<d 


x>0, where: ` 


т — mean objective function parameters vector, 
па! ... M — objective function parameters in each of n years, 
A — input-output coefficient matrix, 
d — the level below which income (profit) cannot drop in 
any year. 


III. INTRODUCTION OF RISK INTO INPUT-OUTPUT 
COEFFICIENTS MATRIX 


Fluctuations of planning parameters are caused either by price 
fluctuation or yield change. Price fluctuations, in terms of LP, affect 
objective function parameters only. On the other hand, yield fluctua- 
tions affect also input-output coefficients. If farm planning is to 
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be consistent with reality, this problem must be taken into account as 
well. 

One of the methods of considering fluctuations of input-output coeffi- 
cients is the so-called ”Chance-constrained Programming” (5). The 
assumptions of this method are the following: if in some constraints there 
are parameters subject to random fluctuations, these constraints cannot 
be met with a 100% probability. To put it in another way, we can assume 
that the risk-affected constraint should be met with a probability of no 
less than for instance 0.90, 0.95 or 0.99. Using the latter approach as the 
starting point, it is necessary to add the 90%o, 950/0 or 99°/o confidence 
interval to the sum of the products of parameters mean values by the 
value of their corresponding activities. Thus, if the deterministic formu- 
lation of the problem is the following: 


a,x > by where: 


ax — vector of input-output coefficients vector, 

x — activity vector, 

b, — the minimum value of ayx ensuring the coherence of the pro- 
gram, 


then it is necessary to replace a,x by: 


t 
e арх — —yVa,x, where 
a 
t — standardized confidence interval, 
Vayx — a,x variation; 


if the constraint is to be met with the required probability. Further: 
Vayx = xTG,x, where: 
б, — variance-covariance matrix of the a, vector. 

The whole equation can thus be presented as follows: 


t — 
a,X— — |/XGyx , 
x y 


and it has to be more than or equal to bx. Then the whole problem is 
as follows: 
maximize: патх 5, such that: 

Ax <b 

ayX—t xTG,x >b 

х z 0. 


5 The objective function has been formulated in a deterministic way to simplify 
the notation, There is no obstacle to formulating it in any other way. 
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Two difficulties are connected with this problem. First, in order to be 
useful, it has to be resolvable and there must be an algorithm of the 
solution. Although this algorithm is available (32), it has a number of 
defects. Not least is its small effectiveness and very high restrictions on 
the size of the problem. The other difficulty with the chance-constrained 
programming is the assumption of a normal distribution of fluctuations 
of input-output coefficients, which is not always tenable. This incon- 
venience can be avoided by using the Tshebyshev inequality (24), in this 
case a considerable increase of the t a parameter has to be taken into 
account.® 

To avoid all these inconveniences, attempts were made to simplify 
this method. Merill (17) and Chen (6) developed methods consisting in 
the interchange of the objective function and the constraint affected 
by input-output parameter fluctuations, when only one constraint is 
subject to tnem. Rahman and Bender (24) formulated a method appli- 
cable in a situation where covariance between input-output parameters 
does not exist or can be ignored. A more general and simplified method 
was worked out by ‘Wicks and Guise (30). It permits the use of LP 
because it is founded on absolute deviation rather than standard deviation 
as a measure of fluctuation. 

Madansky’s method (13) has an entirely different background as it 
is derived from the theory of games. The method assumes that, if there 
is any parameter affected by fluctuations in the constraint, the con- 
straint has to be met in each situation. In terms of the LP used in farm 
organization planning, this means that the constraints under conside- 
ration have to be met each year which is an information source. Thus, 
if the constraint has a deterministic formula: | 


ax < bi, 


then in the case of the a; vector fluctuations and using the Madansky 
assumption, this notation should be presented as follows: 


2 1 
8 The Tshebyshev inequality is: [(Px,—x) < 102] > 1—— which means that the 
t 
probability that the n-th realization will not not deviate from mean by no more 


f е 1 : 
than t times of 6 is higher than 1--- .Thus, if the constraint is to be met with 


1 ZM 
the probability of not less than 1-a, then 1-0 == hence ey —. For 
t a 


the probability equal 0.95, t 24.4, which is more than twice of toos. The smaller 
the a parameter, the bigger that disproportion is. 
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ах < b, where: 
aj — realization of the a; vector in each source-of-information year. 


If b; were also subject to fluctuations, the above notation could be mo- 
dified as: á 


aux < by 


ainX < bin, where: 
by — realization of the b; parameter. 


Introduction of what has been previously achieved into the LP model 
is already obvious: | 4 


maximize: mTx7, such that: 


? As in note 5, 
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IV. SUMMARY: EVALUATION OF METHOD USEFULNESS 


Most of the presented methods have their practical application. The 
methods developed by Freund and by Markowitz are most frequently 
employed (2, 3, 4, 8, 25, 31), but those founded on the theory of games 
are also applied. However, there are no studies whatsoever that would 
compare all the methods in question. More often, we can encounter cri- 
ticism of a particular method, with its weak points and defects 
emphasized. 

A frequent object of criticism is the Freund method. According to 
Petit and Boussard (21), the fundamental objection to the method is that 
it requires an assumption of the normal distribution of yields and prices 
in order to obtain the objective function. This has not been proved so 
far whereas only absolute certainty would justify this assumption. 
Furthermore, there are reports that the distribution of crop yields and 
prices of farm products is not normal or even not symmetric. Petit and 
Boussard after Day (9). Another objection concerns the risk aversion 
coefficient, which is different for every decision-making farmer and has 
to be determined before optimization procedures. This must be determi- 
ned by experiment, which is criticized by Moscardi and de Janvry (18) 
because the coefficient value so defined will be affected by the farmer’s 
attitude towards gambling. 

Although free from the foregoing objections, the Markowitz method 
has also its own defects, the most serious being that a dual solutions is 
not possible (20). 

The above disadvantages of the two methods can be further strength- 
ened by the fact that they require quadratic programming, which is 
more restrictive as to the size of the problem, while the information on 
which the methods are based, that is mean values and variance, is rarely 
credible. In order to obtain such figures, the data covering far more than 
ten years should be used. These data are not always available; moreover, 
the picture can be distorted by yield changes over a longer period due 
to new developments in technology, unless we have the data obtained 
from experiments. Prices can also be affected by such systematic 
changes. 

All these defects of the two methods also hold for their modifica- 
tions, except for that resulting from the use of quadratic programming. 

Methods founded on the theory of games have their own defects as 
well. For example, Wickas and Guise (30) raise an objection that appli- 
cation of the theory of games increases the matrix size. This is an essen- 
tial objection since the LP matrices employed in optimizing farm pro- 
duction and investment plans already have considerable sizes. Another 
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objection Wicks and Guise discuss is that while applying the theory of 
games to risk consideration, we implicitly assume that the farmer’s atti- 
tude towards risk can be described by that theory. There is no evidence 
to support this assumption. Still one more objection can be added that 
information drawn from the past, especially like that used in the theory 
of games, contributes very little to planning. Moreover, with the selec- 
tion of data from previous years, an unconscious assumption is made that 
only those years and none other are representative and their number is 
sufficient as the information basis. 

None of the discussed methods seems to be free from defects. Such 
a method is difficult to imagine, especially until the harmfulness of 
risk is defined. An attempt to deal with the problem in that way is part 
of the safety-first method but it is difficult to apply it in the case of 
fluctuations of input-output coeffients. 


Tab. 1. Crop yields of four oats varieties in lbs. per acre in 1953—1957 


Year 

Variety 1953 , 1954 1955 1956 1957 
A 1472 2112 1920 3520 3072 
B 1568 1984 1824 3104 3328 
C 1440 2368 2496 3552 2848 
D 1952 2688 2184 0t 3200 


ка 


The D crop was destroyed by hail іп 1956. This is the slowest- 
~growing of the four varieties tested. Hail, which normally occurs 
after harvest, affected this variety in 1956 due to a prolonged veg- 
etation period. 


Source: O.L, Walker et al., Application of Game Theoretic Model 
to Decision Making, Agronomy. Journal, no 2, 1964, 
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STRESZCZENIE 


Celem przedstawionej tu pracy jest opis i porównanie metod uwzgledniania 
ryzyka w modelach optymalizacyjnych stosowanych w rolnictwie. Przedmiotem 
opisu były przede wszystkim zagadnienia teoretyczne, a więc zarówno strona for- 
malno-matematyczna prezentowanych metod, jak i formujący je zestaw założeń 
ekonomicznych. i 

Opisywane w niniejszej pracy metody uwzględniania wahań losowych para- 
metrów funkcji celu oparte są na teorii użyteczności bądź na teorii gier, a służące 
do uwzgledniania wahań parametrów techniczno-ekonomicznych również mają uza- 
sadnienie teoretyczne w teorii gier oraz na tak zwanych ograniczeniach losowych 
(chance constraints). Próba oceny wykazała, iż większe nadzieje należy wiązać 
z grupą metod opartych na teorii gier. Do czasu ustalenia na czym polega szkodli- 
wość ekonomiczna ryzyka trudno jest jednak wydawać jednoznaczne oceny. 


РЕЗЮМЕ 


Цель настоящей работы — описать и сопоставить методы учета риска в 
оптимализационных моделях, применяемых в сельском хозяйстве. Предметом 
описания были прежде всего теоретические вопросы, в том чусле как формаль- 
но-математическая сторона представляемых методов, так и формирующий их 
комплекс экономических предпосылок. 

Описанные в настоящем исследовании методы учета случайных колеба- 
ний параметров функции цели опираются на теорию полезности или на теорию 
игр, а методы, служащие для учета колебаний технико-экономических пара- 
метров, теоретически обоснованы также теорией.игр и, кроме того, так назы- 
ваемыми случайными ограничениями (chance constraints). Попытка оценки обна- 
ружила, что больше надежд подае ттруппа методов, опирающихся на сеорию игр. 
Однако до установления, в чем состоит экономическая вредность риска, фор- 
мировать однозначные оценки представляется затруднительным. 
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